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In this work, we introduce an extra singlet pseudoscalar into the Type-III two Higgs doublet model 
(2HDM) which is supposed to solve a series of problems in the modern particle-cosmology. With 
existence of a light pseudoscalar, the h —>■ /rr excess measured at CMS and as well as the {g — 2)^ 
anomaly could be simultaneously explained within certain parameter spaces that can also tolerate 
the data on the flavor-violating processes r —^ /ry and Higgs decay gained at LHC. Within the 
same parameter spaces, the DM relic abundance is well accounted. Moreover, the recently observed 
Galactic Center gamma ray excess(GCE) is proposed to realize through dark matter(DM) pair 
annihilations, and in this work, the scenario of the annihilation being mediated by the pseudoscalar 
is also addressed. 


I. INTRODUCTION 

Even though in the framework of a minimal ex¬ 
tended Standard model (SM) with non-zero neutrino 
masses the leptonic flavor violation (LEV) process is 
negligible for the smallness of neutrino masses which 
are experimentally confirmed [T]. Therefore, a di¬ 
rect search for the LEV processes would provide an 
ideal probe for new physics beyond SM, or in other 
words any observational anomaly may hint us its ex¬ 
istence. Besides the B-factory, since its high energy 
and luminosity, LHC is definitely the machinery for 
the exploration. A search for LEV has been per¬ 
formed by the CMS collaboration via two channels 
h —>■ p,Te and h —>■ and a 2.4 ct excess of the 

branching fraction BR(/i — fir) = (0.84/q; 37)%)[2] 
is reported. If one could account such an excess as 
an anomaly, there should be some mechanisms which 
are obviously beyond SM, to be responsible for it. 
The Type-III two Higgs doublet model (2HDM) is 
one of them, because in the model a flavor-violating 
Yukawa interaction exists which may contribute to 
the LEV at tree level. The model has been explored 
extensively [SE] to study this phenomenological 
observation. Furthermore, the Yukawa interaction 
contributes to the muon g — 2 via one loop dia¬ 
grams and thus would provide a possibility to ex¬ 
plain the {g — 2)^ discrepancy O [HI Il2]. Mean¬ 
while, the flavor-changing Yukawa interaction would 
induce a substantial contribution to the radiative de¬ 
cay T —>■ /ry, thus the flavor-changing Yukawa inter¬ 
action might be rigorously constrained by the avail¬ 
able experimental data isnnuH]. 

One of the main characteristics of the 21st cen¬ 


tury is that the cosmology has already become an 
accurate science and the corresponding observation 
must be combined with the precise measurements 
and new discoveries at the facilities on the Earth to 
testify the standing theories. The identity of dark 
matter (DM) and the interaction which determines 
the behavior of DM particles are the key point and 
searching for them is the most challenging job for 
both experimentalists and theorists of high energy 
physics and cosmologists. 

Recently, the Fermi Large Area Telescope data 
show an excess of gamma-ray at energy of a few GeV 
coming from Galactic Genter (GGE) [T5IET| . To ex¬ 
plain the observation, it is suggested that annihila¬ 
tion of DM particles weighing 30— 70 GeV into bb is 
responsible for the GGE [22H24] . Even though there 
exist other proposals to explain the excess, such as 
a population of millisecond pulsars (MSP) [Ul US] 
which might be responsible for GGE, it is not easy 
to explain the energy spectrum and spatial distri¬ 
bution of the GGE ^7\ HHj . Thus in this work we 
discard the astrophysical source explanation p5H5n] 
and focus on the DM scenario. 

The dwarf galaxies are considered to be the clean¬ 
est sources for detecting gamma rays produced by 
DM annihilations, thus the data on gamma ray ob¬ 
served at Reticulum H pTHHn] may imply existence 
of abundant DM at our galaxy. To be sure, we need 
to compute the cross section of DM annihilation in a 
model of particle physics. Meanwhile other cosmo¬ 
logical phenomena must also be concerned, namely 
the DM annihilation cross section required by the 
new data should be of the same order as that deter¬ 
mined by the thermal DM relic abundance. 

However, the original Type-III 2HDM does not 
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provide a natural explanation to the DM particle 
annihilation cross section. Thus we need to extend 
the model which can also accommodate the DM 
annihilation. It is note that a pseudoscalar could 
mediate the annihilation of dark matter(DM) pair, 
meanwhile due to the small momentum transfer at 
t-channel, the interaction between the DM particles 
coming from the outer space and the nuclei in the 
detector is not affected by the existence of the new 
pseudoscalar, so that the DM particles may evade 
the direct search at not-much sensitive detectors. 
This idea has been implemented in various mod¬ 
els PIISHIH]- It further motivates us to consider 
the DM explanation to the GCE, and introduce a 
Dirac fermion field serving as the DM candidate. In 
this work, we introduce a pseudoscalar oq into the 
Type-III 2HDM. The pseudoscalar does not directly 
couple to the SM particles, but it slightly mixes with 
the CP-odd Higgs which exists in the original Type- 
III 2HDM, thus it would effectively couple to SM via 
this mixing. Therefore the pseudoscalar can mediate 
an effective interaction between the DM xx the 
SM fermions bb. 


Moreover, its introduction may bring up two more 
advantages as following: 

1. An extra pseudoscalar would open a new decay 
channel for the Higgs and thus affect the h —>■ /rr 
excess. 


2. There could be a mixing of the newly intro¬ 
duced pseudoscalar with the CP-odd scalar Aq in 
the original Type-III 2HDM, and hence induces a 
new contribution to the value of (g — 2)^. With in¬ 
crease of the mixing, the contribution of the new 
pseudoscalar would cancel that of the Heavy Higgs, 
thus in this extended model the theoretical predic¬ 
tion on [g — 2)^ can be decreased to a tolerable level. 

In this work the pseudoscalar also plays a role to 
explain the /i —>■ /rr excess and the discrepancy be¬ 
tween theoretical prediction and data of {g — 2)^. 


The paper is organized as follows: after this intro¬ 
duction, we discuss the new scenario where a light 
pseudoscalar is introduced to extend the Type HI 
2HDM in section. [H] In the section. m we investi¬ 
gate the relevant topics including the /i —>■ /ir excess 
observed at CMS, muon {g — 2)^ anomaly, the galac¬ 
tic center gamma ray excess (CCE) and the con¬ 
straints from the t —>■ fij process and LHC Higgs 
data. The numerical results are presented in sec¬ 
tion. IIVI then the last section is devoted to our con- 


II. THE MODEL 

In this work, a Dirac fermion (y) of mass 
stands as the DM candidate and a gauge singlet 
pseudoscalar oq is introduced to extend the Type- 
III two Higgs doublet model, where oq mediates the 
coupling between the dark matter and the SM par¬ 
ticles. The interaction Lagrangian reads 

Aark = -Vx^oXil^X ■ (1) 

The pseudoscalar uq mixes with the pseudoscalar in 
the original 2HDM and then couples to the SM par¬ 
ticles through the potential, given as [15] 

V = V 2 HDM + + portal, (2) 

and 


Vportai = -iBaoH\H 2 + h.c. , (3) 


where i? is a parameter of mass dimension, and the 
Higgs potential is [49] 


H2HDM — + g,2H\H2 -\- {g,^H\H2 -\- h.c.) 

+ Ai (iLjEIi)" + A2 {HIH2)' + A3 {h\h^ 
X {HIH 2 ) - A 4 {h\H 2 ) 


{\^H\H2 + XqH\Hi -\- \'jH\H2) 


(h\H2) 


-I- h.c. 


( 4 ) 


We can explicitly rewrite Hi and H 2 in the Higgs 
basis as 


Hi = 


G+ 

V2 


H 2 = 


H+ 

4'2+iAo 

V2 


( 5 ) 


where G~^ and are the Nambu-Coldstone bosons 
and H~^ and Aq are a charged Higgs boson and a 
CP-odd Higgs boson, respectively. Without losing 
generality, let us concentrate on the CP-conserving 
case, where oq does not develop a vacuum expecta¬ 
tion value (VEV) and all A^ and are set to be real. 
Then the potential is minimized to (Hi) = vlx/2, 
{H 2 ) = 0, with V = 246 CeV. 

In the basis of {(l)i,(j) 2 ), the mass matrix for the 
CP-even Higgs is A4^ whose elements are, 

■^hii — 2Aiu , 

■MI 22 = (6) 

■^hl2 — -^1.21 = ^6^^- 


elusion and discussion. 
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Diagonalizing the matrix, one obtains the physical 
CP-even states h and H {nih < mn) as 


The effective coupling of DM fermions to the medi¬ 
ator given in Eq. Q is simply expressed as, 


4>i 

4>2 


tan 2 Q! = 


cos a 
— sin a 


sma 
cos a 


A/f2 _ yW2 ■ 


with eigen-masses being 




M 


hll + ^h22 


H 

h 


(7) ^dark =-?/x (cos6'a-hsin6<A)xi7^X- (13) 

The Yukawa interactions in the extended Type-III 
2HDM are 

•^Yukawa = ~QlYR ~ QL^CKM^'^Pd ^^R 

- QlHiyWi, - QlH^pl^u'k 

- LlH.yle],, - (14) 


T 



-ML 2 ) +4(A4L2) 


and we consider the eigenstates h and H as the SM- 
like and heavy Higgs bosons respectively. 

The CP-odd Higgs Hq mixes with oq due to the 
potential Tportai (Eq. § , and the mass matrix in the 
(Aq, oo) basis is 


M\ = 


m\ Bv \ 
Bv ml^ ) ’ 


( 8 ) 


where rn\^ = j2 — . Thus, the 

relation between Aq, oq and mass eigenstates A and 
a is 



with the mixing angle 

1 _i, 2Bv , , , 

" = 2 -< >• (“> 


and the masses squares are 






— + AB'^v'^ 


The parameter B in terms of rUa^A and 9 can be 
expressed as 

B = ^{m\-ml) sin2e. ( 11 ) 

The effective coupling of the CP-even Higgs bosons 
to the SM W is igmwC^wwY'' with Chww = 
sin a, Chww = cos a, and the CP-odd Higgs bosons 
A(a) do not couple to W, i.e., ^^( 0 ) 1 ^ 1 ^ = 0. And 
the Hportai is recast in terms of mass eigenstates and 
mixing angle as 

Pportai = “ '^D [■S4e aA + S 2 ff (A^ - Y)] 

X (sin a h-|-cosa H). (12) 


where Q = (ul, Pckm^l)^, L = (Ymns^^l, cl)^ and 
Hi stands for ia 2 H*. Vckm(14ins) is the Cabbibo- 
Kobayashi-Maskawa (Maki-Nakagawa-Sakata) ma¬ 
trix. The general 3-by-3 complex matrices pj induce 
the Higgs-mediated Flavor Changing Neutral Cur¬ 
rent (FCNC). In the mass eigen-basis of the Higgs 
bosons, the Yukawa interactions are recast as 


•^Yukawa — V tjiij fLi4’fRj ^Li{YyuSlS,P^Y ^Rj 

—Ui{VcKMPdPR — pI^ckmPlT^ dj 
-fh.C., (15) 

with (j) = h, H, A,a, f = u, d, e, u, and 


IjAij — 


Vaij — 


m\ 

—- sin a A,- -I- 

11 

Pf 


V2 

m) 

—- cos aSij — 

Pf 


V2 

f -$^os9, 

(/ 

[ %COS0, (/ = 

1 

(/ 

1 

(/ = 


COS a, 


(16) 


where the couplings yAij and yaij exist in the new 
Feynman rules and are accompanied by 75 . For in¬ 
vestigating {g — 2 )^ excess, we do not need to in¬ 
voke the so-called Cheng-Sher ansatz for pj 
since the corresponding parameter space is highly 
restricted mm- The smallness of the mixing pa¬ 
rameter cos a is favored by the current LHC Higgs 
coupling measurements, and we will study the issue 
in later part. In this scenario, the coupling of the 
SM-like Higgs to fermions yhff approaches to the 
SM one, thus the flavor-violating processes mediated 
by the SM-like Higgs boson are almost completely 
suppressed. 
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III. SEVERAL RELEVANT TOPICS 
WHICH ARE SPECIFICALLY ADDRESSED 


proposed by the authors of Ref. [m EH] , the T pa¬ 
rameter in our model is obtained as, 


A. Constraints on the parameter space of the 
afore model 

At first, we explore possible constraints coming 
from B physics and Electroweak precision test, and 
find that the model is more advantageous over the 
Type-II 2HDM as it may evade those constraints 
in the situation of mu+ ^ uia because the tan /3 
enhancement effect does not exist. 


1. Constraints from B Physics 

A light a can mediate the initial and final states 
of decay Bg —)■ in addition to the SM contri¬ 

bution, hence imposes a stringent constraint on the 
model. For rria rnz, the correction due to an a 
exchange at the s-channel was calculated m and 
the results are 


T = 


leTTs^m^ 


F{m\, m^+) cos^ -I- sin ( 


X (E(m^+,m|^) — cos^ , (19) 


with 


^ + y xy X 

F{x,y) = ^ - log-. 20) 

2 x-y y 

The a/A and h/H mixing is highly constrained 
by the current LHC data to be around sin a ~ 1 
and cos6l ~ 1, along with the parameter choice 
mH+ rriA, the T parameter is suppressed, as seen 
in Eq. therefore it does not actually affect appli¬ 
cations of this model as indicated in m- 


B. Relevant processes under investigation 


BR (Rg —>■ yAy ) ~ BR (Rg —>■ y )gj^ (17) 


1 -f 


f{xt,yt,r) 


4m^(m| - rriA) Y (xt) 


for pI^ = 0, and 

BR (Rg y~''y~) « BR (Rg /x+/i“)gj^ (18) 

f{xt,yt,r) 


1-f 


m,Am% - mt 


Y (xt) 


In the following parts, we investigate several rel¬ 
evant processes in this extended Type-III 2HDM 
model. All the puzzles about the h ^ yr excess, 
muon g — 2 discrepancy, constraints coming from 
r —>■ y"f, dark matter relic abundance and GCE ex¬ 
planations which were not solved in previous Type- 
III 2HDM will be addressed. 


1. h ^ yr excess 


for with Xt = yt = 

r = and the / and Y functions can be 

found in Eq. |A.3| of the Appendix. 

The average of the LHCb and CMS measurements 
on this mode is BR(Rg —>• y^y~) = (2.9 ± 0.7) x 
10“® [52U54| . This could be compared against the 
SM prediction, which is taken to be (3.65 ± 0.23) x 
10“® EHlESj- We note that the tan/3 enhancement 
effect in the calculation of Rg —^ y~^y~ with the 
Type-II 2HDM gHl EH [ST] does not appear in our 
model and thus the constraint from Rg leptonic de¬ 
cay is relaxed. 


Existence of the flavor-violating Yukawa coupling 
in the extended Type-III 2HDM may possibly ex¬ 
plain the h ^ yr excess observed by the CMS col¬ 
laboration. Now let us compute the branching ratio 
oih ^ yr in terms of our model, the result is shown 
as 

BR{h ry) = {\yhTn? + \yhnr?) , (21) 

where T/j is the total decay width of the SM-like 
Higgs boson. To meet the observed excess, the flavor 
mixing should be of a magnitude 


2. The T parameter 


p^^ = 



0.0018 



In the Type-III 2HDM, the constraint on the T 
parameter imposed by the electroweak precision test 
may be the most stringent. Following the method 


For T/i ~ 4.2 MeV and cpa = 0.01, as long as 
^ 0(0.1) is reached, the new model is able to 
accommodate the h ^ yr excess. 
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2. muon {g — 2)^ anomaly and t ^ gy 

The previous study indicated that new physics 
may contribute to the anomalous magnetic moment 
of the muon and radiative process r —> gy via a 
chirality flipping dipole operator m 


Crj 

A 2 

^'-NP 


+ h.c., 


(23) 


where i,j denote the flavors of the external lep¬ 
tons, Fap is the electromagnetic strength tensor 
and cr“^ = I [ 7 “,7^]- The diagonal component 
contributes to the anomalous magnetic moment of 
muon, whereas the off-diagonal component corre¬ 
sponds to a dipole transition from t to g. It is 
noted that the coefficient -§r— is derived in var- 

‘^NP 

ious new physics models and has different values 
which would help to determine the corresponding 
parameter space. In this work, we are going to de¬ 
rive this coefficient in the extended Type-Ill 2HDM. 
The flavor-violating Yukawa couplings and the newly 
introduced pseudoscalar induce additional contribu¬ 
tions to {g — 2)^ and r —> gy via one-loop and two 
loop Barr-Zee diagrams m, as shown in the two 
panels of Figure. We include these extra contri¬ 
butions into the numerical computations. 

1. The muon g-2 anomaly 

In our case, the model-dependent coefficient 
C'-^/A^p could be expressed as (SU] 




4m,, 


A 2 

^'-NP 


V 

71' 


(24) 


The contributions coming from the two CP- 
odd Higgs bosons give rise to 




In 


(25) 


Saf=-^^ 




7/ 


(26) 


where (j> = A,a, = mj/m^ and ypff> is de¬ 
fined in Eq.(16), could be found in Eq. A.l 
of the attachm Appendix. The W and Gold- 
stone loops would not contribute to the Barr- 
Zee diagrams for CA(^a)ww = 0- 

2. The T ^ gy process 


With the lepton flavor violating Lagrangian 

£ = emiAf,Ij[ia^^q^{A''lPL + 71)5Pr)/j] (27) 

+h.c, 


the coefficients C*-’ in Eq. (231 could be ex¬ 


pressed in terms of the form factors Ap and 
Ar as done in Ref. [60], 


e Tn-rA^pAj^ 

\/2v 


e rrirAj^pA^r^ 

V2v 


(28) 

Thus, the branching ratio of r —> gy is calcu¬ 
lated as 

BR(r —>■ gy) = BR(t —>■ gvv)x (29) 


t/g 



FIG. 1: The Feynman diagrams for r —>■ gy and {g — 2)^, 
where (j) represents h,H, A, a and the loop at the bottom 
panel corresponds to t, W and Goldstone loops. 

Due to lyfiTul — Ivtiprl in the model, we have 
\Al\ = \Ar\ = \A\. The two CP-odd Higgs 
bosons contribute to the form factor A through 
1-loop and 2-loop Barr-Zee diagrams as shown 
in Fig. 12 and the new form factor is obtained 
as 

1 


A = 


IbTT^ 


(-^1 + -^ 7 ) • 


( 30 ) 


where 


Ai = V 2 J 2 


y4>^Ty4>T 






At,b _ NcQ}o^em 

* * TT - 


p-f 
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The other contributions from the CP-even 
Higgs bosons are induced by the 2-loop Barr- 
Zee diagrams where an intermediate photon 
and a bP-boson are involved [THUS]- 


3. DM annihilation and GCE 

For TOa ^ mA, the Dirac DM fermions annihilate 
into bb primarily through exchanging a at s-channel, 
the annihilation cross section for the relative velocity 
Vr is given as 


(avr) ~ 


yliply 


’20' 


IGtt (s — -I- ' 


(33) 


where y/s is the center-of-mass energy of the anni¬ 
hilating DM fermions, rua and Fq are the mass and 
decay width of the mediator boson a respectively. In 
the non-relativistic approximation s 


thus Eq. 33 can be rewritten as 




r)^ 


vliPj 


^29 


2567rm2 (J -|- -p 


(34) 


where 7 = iriaTa/'im'^ and 6 are two dimensionless 
parameters, and the kinematic factor S is defined as 
(5 = 1 — m^/(4TO^). 

As long as i5 is not too small, the DM annihi¬ 
lation could occur in the region far away from the 
resonance, then the cross section is almost indepen¬ 
dent of the velocity. In this case, the GCE and cor¬ 
rect thermal DM relic density could be accommo¬ 
dated simultaneously provided that the parameter 
8^0 is adjusted to an appropriate value. Eor 
small (5, the resonance effect would enhance the DM 
annihilation cross section. Eor that case, adjusting 
the parameter y\{p'‘f )‘^S 20 /{5'^ + 7 ^) can give a rea¬ 
sonable explanation to the GCE observation. When 
(5 > 0, Eq. 34 indicates that the magnitude of {aVr) 
decreases as the temperature increases and the pro¬ 
cess XX —>■ a —does not sufficiently reduce the 
DM abundance at the freeze-out epoch, therefore 
some other annihilation processes which affect the 
DM relic abundance must exist in the Lee-Weinberg 
evolution equation. 


IV. A SYNTHESIS OF ALL THE 
INGREDIENTS 

In this section, we perform a numerical analy¬ 
sis to investigate the CMS /i —)■ /rr excess, muon 


{g — 2)^ anomaly and the r —>■ /rq, as well as DM 
relic abundance and GCE in the extended Type- 
III 2HDM. The model has been implemented in the 
program FevnRules[63j. and the model file of the 
form CalcHEP [M] has been employed in the pack¬ 
ages micrOMEGAs 4.1.8 [SS] to calculate the relic 
density and the annihilation cross section of DM. 

The Higgs masses in the model are set as: mu = 
150 GeV, mH± = ruA = 300 GeV. Here the value of 
mH± is allowed by the flavor physics constraint [ 66 ] , 
and niA = itih± is suggested by the T parameter 
constraint as indicated in section. HI A 2 and the 
value of niH employed in our computations is con¬ 
sistent with the vacuum stability requirement [67]. 
The magnitude of the two CP-even Higgs mixing 
angle a, is severely constrained by the recent Higgs 
data [ 68 ], and its closeness to 7r/2 will be explored 
in later parts of this work. 


The invisible and undetected decays of SM-like 
Higgs boson are accounted as the decays of beyond 
SM (BSM) Higgs boson [SH], and it is denoted as 
r(/i —>• aa) + r{h —>• gr) = FBSM(/i)- Since the 
contribution of BSM to the decay branching ratio is 
bounded bellow 0.34 at 68 % CL, the decay oih ^ aa 
and h —> /rr in this model would be constrained. 
The best fit of the branching ratio of the undetected 
decays is < 0.23 at 68 % CL, and this limit allows 
that of BR{h —)• gr) = (0.84^q 37 )%, therefore it 
implies that the constraint on the flavor violating 
process is relaxed. As ma < mh/ 2 , the a/A mixing 
angle 9 dominates the exotic decay rate of the SM- 
like Higgs h ^ aa and changes the total decay width 
of the SM-like Higgs, by which the prediction value 
of BR(/i —>■ p,r) in this model would be affected. 
The magnitude of 9 is required to be at order of 
0(0.01 — 0.1) required by the present Higgs signal 
fit which is also welcome by the estimate of the DM 
relic abundance and GCE interpretation since this 
value affects the magnitude of {avr){ see Eq|3|). 

The parameter dependence of the coupling 
denoted as pbt for simplicity) needs to be carefully 
analyzed for the following reasons: the coupling 
Pbb is responsible for one of the dominant decay 
modes of the SM-like Higgs boson h —>■ bb, thus 
any change of pbb would affect the theoretical pre¬ 
diction on BR(/i —>■ pr); meanwhile the DM annihi¬ 
lation XX bb is supposed to be the dominant one 
and its cross section is proportional to the square of 
Pbb, as given in Eq[^ pbb is set at the same order 
as Pf^r, i-e. as aforementioned in secjlll B 1| The 
contributions of tau lepton and top quark dominate 
BR(r —)• /rq), as shown in Eq|M] and Eq|^ Thus 
the two parameters pl'^ and (denoted as prr and 
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Pit for simplicity) need to be explored, here we use 
the parameter range of Ptt^Pu) ~ Pbb- Even though 
we do not take the Cheng-Sher ansatz for , we still 
choose a negative value for as in the Cheng-Sher 
ansatz lana, except for prr and pu while consider¬ 
ing the current experimental constraint of r —>■ py. 

To interpret GCE, the DM fermion mass and the 
coupling of the DM fermion to the mediator (a, A) 
are set as rrix = 30 GeV and = 0.5, and the range 
of rua ranges between 30 GeV-95 GeV to account for 
the resonance effects in the DM annihilation process. 

To obtain the parameter spaces favored by the 
physical picture including all the aforementioned 
constraints, we perform a complete numerical anal¬ 
ysis for all possible parameter spaces: a — p^r, a — 0, 
9 - Pf^r, Pbb - PfiT, Pit - Ptt, and nia - 9, with rele¬ 
vant parameters being free within the ranges of: 30 
GeV< ma < 95 GeV, 0.025 < 9 < 0.1, 1.475 < 
a < 1.57, —0.115 < Pfj,r < 0, —0.3 < pbb < 0, 
—0.28 < Pu < 0.28, —0.05 < prr < 0.05 based on 
aforementioned arguments. For each specific param¬ 
eter space listed above, several parameters need to 
be fixed as shown in the following 

• for parameters spaces of a — 9 (a — p^r)'- 

Pp.1^ = - 0 . 01 , Prr = 0 . 012 , Ptt = - 0 . 2 , Pbb = 

-0.2, and ma = 50GeV(46 GeV), p^^ = 
-0.102(6» = 0.06) ; 

• for parameters spaces of pu — Ptt (pbb — PtiT)- 
ma = 50GeV, a = 1.546, p^^ = —0.01, 9 = 
0.045, Pbb = —0.2 and p^^ = —0.102 {p^-T = 
0.012 and pu = —0.2); 

• for the PfiT — 9(ma — 9) parameter spaces: 

a = 1.546, Ptt = 0.015, pu = -0.2, pbb = 
—0.2, and p^^ = —0.01(—0.02), TOq = 50GeV 

{PtaT = - 0 . 102 ). 

Conducting a numerical analysis by means of the 
above parameter setup, the relevant processes are 
depicted by the figure (see the caption of the figure 
for details). 


1. The h ^ p,T 

To understand the CMS excess BR(/i —>■ pr) = 
0.84 +q;| 7%, the relation between the coupling 
Pt^T and a given in Eq. 22 should be satisfied. 
The region colored by magenta in the top-left 
plot of Fig. is allowed to explain this ex¬ 
cess. With a properly fixed a, the both plots of 
the middle panel show that a sizeable coupling 
PtiT ~ 0.1 is required. 


The exotic decay mode h ^ aa with a large 
branching ratio can efficiently change the to¬ 
tal Higgs decay width, thus the ajA mixing 
angle 9 and the mediator mass ma affect the 
prediction of BR(/i —)• pr), as plotted in the 
bottom-right panel. 

The whole magenta region in the bottom-right 
panel with the value of BR(/i —>■ pr) = 0.61 
implies that when 9 is relatively small or in a 
case of 2 toj(, < m^, the decay process oih ^ aa 
cannot occur, so would not affect the Higgs 
total width and becomes irrelevant to —> pr 
process. Eqj^ clearly interprets the situation. 

It is worth indicating that when this work is 
close to be finalized, the ATLAS collaboration 
published new analysis on —>■ pr^,, they ob¬ 
tained a slightly smaller excess compared to 
the CMS result, while its upper bound is con¬ 
sistent with the CMS result [55] ■ 

2. The muon g — 2 anomaly 

The contribution of two-loop Barr-Zee dia¬ 
grams to SOfa is negligible because of the small¬ 
ness of cos a. The dominant contributions to 
SOfa include the one-loop diagrams where CP- 
even Higgs H and CP-odd A are mediators, 
especially, the two diagrams respectively pro¬ 
vide negative and positive contributions. 

The contribution of the one-loop diagram 
where CP-odd a is the mediator to Sa^i is pos¬ 
itive, it becomes larger for smaller ma and/or 
larger 9, therefore cancels out the contribution 
of the diagram where H is involved. The situ¬ 
ation is depicted at the top- and bottom-right 
panels of Fig. (by the red dashed curves). 

To gain the outcome result which is con¬ 
sistent with the present experimental value 
(5a^/10“® = (2.61 ± 0 . 8 ) (ZOI), a sizeable 
is required, it is also consistent with the CMS 
h ^ pr excess, and its dependence on p^^ 
is demonstrated in the plots(top-left and the 
middle panels). 

Within the allowed parameter space in p^-T- 
Ptt, a value of muon g-2: 8a^ = 2.80 x 10 “® is 
reached. 

3. Constraints from t —>■ py 

The dominant contributions come from the 
one-loop diagrams (A-loop and H-loop di¬ 
agrams) and the Barr-Zee diagram with 
top quark being involved, as shown in 
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0.03 
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FIG. 2: The corresponding results are presented as functions of a, pbb, ptt, Ptt, P/it, ma and 6. The cyan (orange) 
regions give the DM annihilation cross section of (aVr) = (0.5 — 4) x 10“^®cm®/s((0.5 — 1) x 10“^®cni®/s). DM relic 
density Q,h?{— 0.1197±0.0022) is depicted by blue lines and contours. Magenta region/contours are for CMS h —>■ pr 
excess. Red dashed lines are contours of 5a^/10“®. The black curves are contours of BR(r —>■ p'y)/lQ~^, and the red 
region in the bottom-left panel is consistent with the experimental constraint < 4.4 x 10“®. The green curves are 
contours of BR(/i —>■ aa). The gray region is ruled out by the constraints of CMS Higgs signal strengths. 


Eq. ( |30|31|32 |. Existence of opposite signs be¬ 


tween the contributions of the CP-even Higgs 
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H (negative) and the CP-odd Higgs A (posi¬ 
tive) one-loop diagrams [T^] leads to a cancela¬ 
tion effect, however, it does not occur as long 
as tua > ran- With a large coupling pu and 
under the limit sin a ~ 1, the Barr-Zee dia¬ 
gram involving top-quark composes the domi¬ 
nant contribution (positive) to r —>■ pj. 


Additionally, like the case of i5a^, as rria is 
smaller and 9 is larger, the importance of the 
one-loop where a is involved to BR(t —>■ pj) 
enhances as shown in the bottom-right panel 
of Fig. As well, a smaller a would also in¬ 
duce an enhancement of the contribution of the 
h mediated Barr-Zee diagram with a W loop. 
Thus due to the contributions of this Barr-Zee 
diagram and the a one-loop diagram, for a re¬ 
gion with larger 9 and smaller a, our calcula¬ 
tion would predict an even larger BR(t — pj) 
(see the top-right panel). 

The theoretical prediction on BR(t — p'y) 
with respect to Pi^r (plotted at the top-left and 
middle panels) shows that its behavior is sim¬ 
ilar to that for (5a^. 


It is noteworthy that since the the form fac¬ 
tor A is related to tau lepton and to p qua rk 
through the concerned loops, (see Eq. ( |31|32 )), 
the current experimental bound BR(r 
p-f) < 4.4 X 10"® [ZH [ig constrains prr and 
Pit strictly. The red region in the bottom-left 
panel of Fig. is allowed. It is interesting 
to note that ptt and prr should have opposite 
signs as favored by the data. This result agrees 
with that given by the authors of Ref. [Ill- 
One can see that the upper bound demands 
p^^ to be small as about \prr \ < 0.04. 


values of (avr) = 1.8(2.5) x 10 ^®cm®/s. 

The blue curves and contours from Fig. [^rep- 
resent the correct dark matter relic density 
(0/i2 = 0.1197 ± 0.0022 [ig). We note that 
the range of {aVr) favored by the DM relic den¬ 
sity and GCE as expecte d highl y depends on 
nia, 9, and ptb (see Sec. IIIB3). Due to the 
enhancement of the annihilation cross section 
when the mediator mass is close to ~ 2 to^, 
the dark matter relic density rules out a range 
in the parameter space (see the bottom-right 
panel). 


With a sizeable pbb ~ 0(0.1), the DM relic 
density and GCE could be tolerated simul¬ 
taneously in a range of (aVr) = (0.5 — 4) x 
10 “^®cm®/s. 


However, the results newly reported by the 
Fermi-LAT and DES collaborations constrain 
(aVr) to be smaller than 1.0 x 10“^®cm®/s 
for the 30 GeV dark matter fermions, which 
is slightly smaller than the value required by 
the thermal relic abundance, thus there should 
exist other additional DM annihilation chan¬ 
nels to make up the correct DM relic den¬ 
sity [751175] . 


5. The CMS constraints 

Recently, the CMS collaboration combined 
the comprehensive sets of production and 
decay measurements for the 125 GeV- 
Higgs boson, including decay channels into 
j'y, ZZ* ,WW* ,T'^T~, bb, and p'^p~, and 
found no significant deviation from the stan¬ 
dard model predictions |5H]. That synthesis 
should severely constrain the parameter spaces 
of all built models. 


4. DM relic density and GCE 

As long as assuming that the observed GCE is 
caused by dark matter annihilation, the lower 
bound of the annihilation cross section {aVr) 
should be about 0.5 x 10“^®cm®/s as discussed 
in Ref.jZg. The upper bound of {aVr) is deter¬ 
mined to be 4.0 X 10“^® cm®/s at 95% CL. [33] . 
Also the data of Pass 8 of Fermi-LAT [73| from 
dwarf spheroidal satellite galaxies set a new 
upper bound (~ 1.0 x 10“^®cm®/s) on the dark 
matter annihilation cross section at the dark 
matter mass 30 GeV. The areas in Fig. [^ that 
give an annihilation cross section of {aVr) = 
0.5 — 4.0(1.0) X 10“^®cm®/s are depicted as 
the cyan (orange) regions. As for the top-left 
(bottom-left) parameter space, the calculated 


More specifically, the signal strengths for 
T+r“, 66, LFLF*, ZZ*, 77 channels are defined 
as 


<j{h) X BRi 
crSM(/j) X BRf“ 


(35) 


which can be gained by fitting the CMS data 
within Icr tolerance [77]. In the text (j(h) and 
(t®'^(6 .) ( BRi and BR™) correspond to the 
Higgs production cross section (decay branch¬ 
ing fractions of the five decay modes) pre¬ 
dicted by respectively this model and the SM 
as BRj = T{h —)• ii)/ (rBSM(M + rsM(^))- 
We find that the upper bound on pw and 
the lower bound on pj set more rigorous con¬ 
straints on the parameters of this model, and 
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CMS[68] 

nia — 9 

a — 0 


PfLT G 

pTT Ptt 

P 1 J.T Pbb 

(JiT 

0.91 ±0.28 

0.94-1.07 

0.92-1.11 

0.93-1.10 

0.94-1.08 

0.85-1.19 

0.92-1.0 

Pb 

0.84 ± 0.44 

0.58-0.67 

0.63-0.65 

0.64-0.65 

0.58-0.67 

0.64-0.69 

0.64-0.66 

pw 

0.83 ±0.21 

0.88-1.0 

0.88-1.04 

0.89-1.03 

0.88-1.0 

0.96-1.04 

0.87-1.03 

pz 

1.0 ±0.29 

0.88-1.0 

0.88-1.04 

0.89-1.03 

0.88-1.0 

0.96-1.04 

0.87-1.03 


1.12 ± 0.24 

0.88-1.0 

0.88-1.04 

0.89-1.04 

0.88-1.01 

0.97-1.04 

0.88-1.04 


TABLE I: A comparison of the signal strengths in the allowed parameter spaces of this model with the Higgs data 
coming from CMS [68] . 


the grey areas in Figj^are excluded. With the 
allowed spaces of the parameters, the ranges 
of the signal strengths Hi are given in Table. |Tj 

The CMS Higgs decay signal strength strongly 
constrains the h/H, a/A mixing angles a and 
9, concretely it demands both sin a and cos 9 to 
be close to unity (see the middle-left and top 
panels). The constraint on a is found to be 
compatible with the parameter space obtained 
by accounting for the h ^ excess, ptb is 
also severely bounded since it could affect the 
primary decay mode of the SM-like Higgs sig¬ 
nificantly. As is shown in middle-right panel 
of Fig. which is a contour diagram of ptb- 
one notes that only a small band with pbb 
being between —0.15 and —0.22 is not ruled 
out. 

Meanwhile from the contour diagram of rua—d, 
one can see that as long as BR(/i —)■ aa) is 
larger than 0.12, the grey regions in the fig¬ 
ure are excluded by the CMS Higgs signal 
strengths. 

For the not-yet-detected channel h ^ aa , 
the pseudoscalar primarily decays into bb, then 
there should be an additional contribution of 
the process h ^ aa ^ Ab in the h ^ bb 
searches [75], and the data of the CMS exper¬ 
iment [Ml [79] would definitely constrain the 
coupling between hbb as long as rua < m;j/2. 
The upper bound on the undetected decay of 
Higgs will be further improved as the bb pair 
production is measured at the 13 (14)TeV. 
Moreover the h —>■ 2b2v searches at LHC |25| 
could give more rigorous constraints on the pa¬ 
rameters nia and 9. 

6 . Collider observation of the pseudoscalar 

For the benchmark scenarios of the work, when 
TTia > 2.m^, the masses relation > rnt+ma 
opens the channel pp{gg) —>■ A —> h{h —>■ 
77)0(0 — )■ xx ) for the observation of the pseu¬ 


doscalar o at the LHC, which is found to be the 
case of the mono-Higgs searches in [SD]. The 
backgrounds are dominated by the SM pro¬ 
cess pp —>■ with Z w, and the Higgs 

associated production process pp —>■ Zh with 
Z —>■ vv. The collider analysis is performed by 
generating signal and background events with 
MadGraph5_aMC@NLO HU Eg at 14 TeV, 
and then passing on to Pythia 8.1 [84] for par- 
ton shower and hadronization, and the detec¬ 
tor simulation is conducted by Delphes 3 [55] 
at last. After implement the selection cuts 
rury^ S [120,130] GeV, as well as PfC > 76 
GeV following Ref. [50] , the number of events of 
signal and background are obtained as 37 and 
48 respectively, thus the significance for the 
observation of a is found to S/y/S + B ^ 4 for 
the benchmark of rUa — 9 with 9 = 0.08, rUa = 
76 GeV and with an integrated luminosity 
C = 300/6-1. 

The probe of the pseudoscalar could also be 
conducted with hard b-jets and large missing 
transverse energy [86] when rUa > 2m^ and the 
pseudoscalar dominantly decay to dark mat¬ 
ter. And when the pseudoscalar decay mostly 
to b-quarks, the collider search of the pseu¬ 
doscalar could be found in Ref. [71] . 


V. CONCLUSIONS AND DISCUSSIONS 

The discovered anomalies by the LHC experi¬ 
ments and observations of the Dark matter greatly 
excite the curiosity of human beings and inspire en¬ 
thusiasm of searching for new physics beyond stan¬ 
dard model. However, so far, the trend is not 


^ Here, we would like to mention that including the back¬ 
ground Z^ + jets (with a jet faked a photon) might reduce 
the observation significance a bit 
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very successful. Even though we know new physics 
must be around, but do not know its scale. Direct 
and indirect search for dark matter, LHC experi¬ 
ments, long-baseline and short-baseline neutrino ex¬ 
periments and numerous lower energy experiments 
including BELLE, BES and many others provide 
hints towards new physics beyond standard model, 
however at the same time, set more and more rig¬ 
orous constraints on those models which have been 
proposed to explain the anomalies observed in as¬ 
tronomy and high energy experiments. Some of the 
models survive the so-far measurements and many 
have been ruled out. 

In this work, we have extended the Type-Bl 
2HDM by introducing a pseudoscalar a which can 
mediate the DM pair annihilation process. In a re¬ 
cent paper, Han et al.[S7] also extended the 2HDM 
with an aligned Yukawa sector to explain the {g—2)^ 
excess, in comparison, our scheme is somewhat dif¬ 
ferent from theirs. We not only consider the {g — 2)^ 
excess, but also many other constraints from both 
earth experiments and cosmology. In this frame¬ 
work, the LEV process h ^ gr observed at LHC 
is addressed, the possibility of explaining the muon 
g—2 anomaly, dark matter relic abundance and GCE 
have also been investigated, the role played by the 
newly introduced a is studied in some details. It 
is found that there indeed exist certain parameter 
spaces which can tolerate h —> /rr excess, muon g — 2 
discrepancy and dark matter relic abundance. 

With the flavor violating coupling ^ 0(0.1) 
and a tiny mixing between h and H around sin a ^ 
1 , the observation of the h ^ pr excess could be 
easily accommodated. The pseudoscalar a opens an 
important undetected decay channel for the SM-like 
Higgs {h —>• aa) as rua < m^/2, thus affects the 
branching ratio of h —> rp. It also plays a role to 
explain the discrepancy between theoretical predic¬ 
tion and data of {g—2)^. A smaller Wq and a slightly 
larger CP-odd a/A mixing 9 helps to interpret the 
muon g — 2 anomaly. Increasing the CP-odd Higgs 
mixing angle and \p^r \ increases BR(t —)■ /ry), thus 
would further constrain the parameters space of our 
model. The measurement of the branching ratio of 
T —^ /ry sets a stringent bound on the flavor con¬ 
serving Yukawa couplings ptt and p^r^ moreover, it 
determines an opposite sign between the two cou¬ 
plings. There are parameter regions in our model 
allowed by the current experimental data of r —>■ /ry 
where both the measured h —>■ pr excess and the 
muon g-2 anomaly can be explained. 

To account for the dark matter relic density, the 
a/A mixing angle should be of order 0(0.01) and 
Pbb 0(0.1) is required, however this parameter 


region does not coincide with that favored by the 
GCE observation if it is postulated that the GCE 
is fully coming from the XX DM pair annihilation. 
This inconsistency may imply that there are other 
sources to result in the GCE besides the pure %% 
DM annihilation mechanism. 

The LEV process h —> /rr is studied in terms of our 
model and our prediction on BR(/i —>■ pr) is quali¬ 
tatively consistent with that observed by CMS and 
ATLAS at 8 TeV, however to make a decisive conclu¬ 
sion, more data are needed and the LHC Run H of 14 
TeV should help. A synthesis of data accumulated 
by high energy collider LHC, the future SPPC of 
100 TeV and maybe some lower energy experiments 
as well as the new astronomical observation would 
make the whole picture clearer, then we will be able 
to judge whether this model indeed works or needs 
to be further modified. For the case of rUa > 2m^, 
the mono-Higgs search provide one possible probe of 
the pseudoscalar. 


Appendix: Loop functions 


The function and h using for the calculation of 
g — 2 and t ^ pj are given by. 


fA,air) = 


1 1 1 2 :(I - x) , . 

dx—log- (A. I) 


(I-x)- 


The functions / (x, y, r) and Inami-Lim function 
Y{x) being used in Eq. 0 are 


fix,y,r) = - 


r (x — 1) — X , x\ogx 

logr ' 


(r-I) (x-I) 
y\ogy X log y 


(y-I) {r-x){x-l) 

, ^ , X a; — 4 , 3a: log x 


(a: - I) 

, (A.2) 

(A.3) 
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